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Abstract 

Recently a lot of new experimental results on open charm hadrons have appeared. In 
particular many D meson resonances have been discovered. We discuss strong decays 
of positive and negative parity charmed mesons within heavy meson chiral pertur- 
(Nj \ bation theory and study the impact of excited charm states on the determination of 

the effective meson couplings pQ. Motivated by recent experimental results we also 
■ reconsider semileptonic D — > Plv\ and D — > Vlv\ decays within a model which com- 

bines heavy quark symmetry and properties of the chiral Lagrangian. Using limits 
of soft collinear effective theory and heavy quark effective theory we parametrize the 
semileptonic form factors. We include excited charm meson states in our Lagrangians 
and determine their impact on the charm meson semileptonic form factors. In some 
scenarios of new physics an up-like heavy quark appears, which induces FCNC at 
tree level for the c — > uZ transitions. We investigate FCNC effects in D rare de- 
cays in particular the c — > ul + l~ transition which might occur in D + — > 7r + Z + Z~ and 

D° -> (Pin-. 
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1 Strong decays of positive and negative parity 
charmed mesons 

The strong and electromagnetic transitions of positive and negative parity charmed 
mesons have already been studied within a variety of approaches (see references [8] - 
[23] given in [TJ). In ref. [2] the chiral loop corrections to the D* — ► Dtt and D* — ► D^j 
decays were calculated and a numerical extraction of the one-loop bare couplings was 
first performed. Since this calculation preceded the discovery of even-parity meson 
states, it did not involve loop contributions containing the even-parity meson states. 
The ratios of the radiative and strong decay widths, and the isospin violating decay 
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Figure 1: Sunrise (left) and sunrise road (right) topology diagrams 

D* — > D S 7T° were used to extract the relevant couplings. However, since that time, 
the experimental situation has improved and therefore we consider the chiral loop 
contributions to the strong decays of the even and odd parity charmed meson states 
using HH^PT. In our calculation we consider the strong decay modes D* + , D*°, Dq + 
D*° and D[° (given in Table 1 of [I]) 

The existing data on the decay widths enable us to constrain the leading order 
parameters: the D*Dix coupling g, DqDh coupling h, and the coupling g which enters 
in the interaction of even parity charmed mesons and the light pseudo-Goldstone 
bosons. Although the coupling g is not yet experimentally constrained, it moderatelly 
affects the decay amplitudes which we investigate. 

Due to the divergences coming from the chiral loops one needs to include the 
appropriate counterterms. Therefore we construct a full operator basis of the relevant 
counterterms and include it into our effective theory Lagrangian. The details of the 
Heavy hadron chiral perturbation theory HH^PT we use is given in pp. First we 
determine wavefunction renormalization of the relevant heavy meson fields considering 
the effects of the chiral loops given by the left diagram in Fig. [TJ Then we calculate 
loop corrections for the PP*rr, PqPiTt and PqPh vertexes. At zeroth order in 1/mq 
expansion these are identical to the P*P*7r, P*P*ir and PlP*n couplings respectively 
due to heavy quark spin symmetry (right diagram in Fig. [T]). 

Using known experimental values for the decay widths of D + *, Dq*, Dq* and D l , 
and the upper bound on the width of D°* one can extract the values for the bare 
couplings g, h and g from a fit to the data. The decay rates are namely given by 

I eff. 12 

r^^^H^^M 3 , (i) 

and a similar expression (up to polarization averaging phase space factors) for T(P ~~ * 
n P) and r(P* — ► 7r P*) with g coupling replaced by h and | fc^ | 3 replaced by E^.\k n \. 
Here k n is the three-momentum vector of the outgoing pion and E n its energy. The 
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Calculation scheme 
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\h\ 
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Leading order 


0.61 


0.52 


-0.15 


One-loop without positive parity states 


0.53 






One-loop with positive parity states 


0.66 


0.47 


-0.06 



Table 1: Summary of our results for the effective couplings as explained in the text. 
The listed best-fit values for the one-loop calculated bare couplings were obtained by 
neglecting counterterms' contributions at the regularization scale fi ~ 1 GeV. 



renormalization condition for the couplings can be written as 



9p*p b ^ ~ 9 , 9^p* Pb ^ (/J 



Z\p a p*^ 



with similar expressions for the h and g couplings. 

We perform a fit with a renormalization scale set to \i ~ 1 GeV [2] and we neglect 
counterterm contributions altogether. Our choice of the renormalization scale in 
dimensional regularization is arbitrary and depends on the renormalization scheme. 
Therefore any quantitative estimate made with such a procedure cannot be considered 
meaningful without also thoroughly investigating counterterm, quark mass and scale 
dependencies. We perform a Monte-Carlo randomized least-squares fit for all the 
three couplings in the prescribed regions pQ using the experimental values for the 
decay rates to compute x 2 an d using values from PDG [6] for the masses of final state 
heavy and light mesons. In the case of excited Dq and D[ mesons, we also assume 
saturation of the measured decay widths with the strong decay channels to ground 
state charmed mesons and pions (Dq — > Dtt and D[ — > D*n). 

Due to the rather large mass splitting between positive in negative parity states 
Ash, we find that the perturbative expansion holds for scales below n < 1 GeV, while 
these new strongly scale dependent corrections become large at higher renormalization 
scales. From the data for the four decay widths we obtain the best-fitted values for 
the bare couplings, which we summarize in Table [TJ We are able to determine all 
the three couplings since the contributions proportional to the coupling g appear 
indirectly, through the loop corrections. 

Since we consider decay modes with the pion in the final state, we do not expect 
sizable contribution of the counterterms. Namely, the counterterms which appear in 
our study are proportional to the light quark masses, and not to the strange quark 
mass [2]. The effects of counterterm contributions in the decay modes we analyze, are 
nevertheless estimated by making the random distribution of the relevant counterterm 
couplings (see [1]). The counterterm contributions of order 0(1) can spread the best 
fitted values of g, \h\ by roughly 15% and g by as much as 60%. Similarly, up to 20% 



3 



shifts in the renormalization scale modify the fitted values for the g and \h\ by less than 
10% while g may even change sign at high renormalization scales. Combined with the 
estimated 20% uncertainty due to discrepancies in the measured excited heavy meson 
masses, we consider these are the dominant sources of error in our determination of 
the couplings. One should keep in mind however that without better experimental 
data and/or lattice QCD inputs, the phenomenology of strong decays of charmed 
mesons presented above ultimately cannot be considered reliable at this stage. 

A full calculation of the strong decay couplings should also contain, in addition 
to the calculated contributions, the relevant l/rrin corrections as discussed in ref [3]. 
There, the next to leading terms (l/m#) were included in the study of charm meson 
mass spectrum. Due to the very large number of unknown couplings the combination 
of 1/rriH and chiral corrections does not seem to be possible for the decay modes we 
consider. In addition, recent studies of the lattice QCD groups [U [5] indicate that 
the \ I van corrections do not contribute significantly to their determined values of the 
strong couplings, and we therefore assume the same to be true in our calculations of 
chiral corrections. 

Due to computational problems associated with the chiral limit, lattice QCD 
studies perform calculations at large light quark masses and then employ a chiral ex- 
trapolation — > of their results. Our analysis of such chiral extrapolation of the 
coupling g shows that the full loop contributions of excited charmed mesons give siz- 
able effects in modifying the slope and curvature in the limit — > 0. We argue that 
this is due to the inclusion of hard pion momentum scales inside chiral loop integrals 
containing the large mass splitting between charmed mesons of opposite parity Ash 
which does not vanish in the chiral limit. If we instead impose physically motivated 
approximations for these contributions - we expand them in terms of 1/Ash - the 
effects reduce mainly to the changes in the determined values of the bare couplings, 
used in the extrapolation, with explicit h contributions shrinking to the order of 5%. 
Consequently one can infer on the good convergence of the 1/ Ash expansion. 

As a summary of our results, we point out that chiral loop corrections in strong 
charm meson decays can be kept under control, but they give important contributions 
and are relevant for the precise extraction of the strong coupling constants g, h and 

g- 

2 Charm meson resonances in D semileptonic de- 
cays 

The knowledge of the form factors which describe the weak heavy — > light semilep- 
tonic transitions is very important for the accurate determination of the CKM param- 
eters from the experimentally measured exclusive decay rates. Usually, the attention 
has been devoted to B decays and the determination of the phase of the V u b CKM 
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matrix element. At the same time in the charm sector, the most accurate determi- 
nation of the size of V cs and V c d matrix elements is not from a direct measurement, 
mainly due to theoretical uncertainties in the calculations of the relevant form factors' 
shapes. 

Recently, there have been new interesting results on D-meson semileptonic decays. 
The CLEO and FOCUS collaborations have studied semileptonic decays D° — > 7r~£ + z/ 
and D° — ► K~i + v [3, E]. Their data provide new information on the D° — ► n~£ + v 
and D° — > K~l + v form factors. Usually in D semileptonic decays a simple pole 
parametrization was used in the past. The results of Refs. [7J E] for the single pole 
parameters required by the fit of their data, however, suggest pole masses, which are 
inconsistent with the physical masses of the lowest lying charm meson resonances. 
In their analyses they also utilized a modified pole fit as suggested in [9] and their 
results indeed suggest the existence of contributions beyond the lowest lying charm 
meson resonances [7j. 

In addition to these results new experimental studies of charm meson resonances 
have provided a lot of new information on the charm sector [HE [TTj [T2l [T3] which we 
can now apply to D and D s semileptonic decays. 

The purpose of our studies [TJ1 [151 IB] is to accommodate contributions of the 
newly discovered and theoretically predicted charm mesons in form factors which are 
parametrized using constraints coming from heavy quark effective theory (HQET) 
limit for the region of g^ ax and in the q 2 ~ region using results of soft collinear 
effective theory (SCET). We restrain our discussion to the leading chiral and l/m# 
terms in the expansion. 

The standard decomposition of the current matrix elements relevant to semilep- 
tonic decays between a heavy pseudoscalar meson state \H{p H )) with momentum p v H 
and a light pseudoscalar meson state \P(pp)) with momentum p p is in terms of two 
scalar functions of the exchanged momentum squared q 2 = (pn — Pp) 2 - the form 
factors F + (q 2 ) and Fo(q 2 ). Here F + denotes the vector form factor and it is dom- 
inated by vector meson resonances, while F denotes the scalar form factor and is 
expected to be dominated by scalar meson resonance exchange [TTKTS]. In order that 
the matrix elements are finite at q 2 = 0, the form factors must also satisfy the relation 
F+(0) = F o (0). 

The transition of \H(pu)) to light vector meson \V(py,ev)) with momentum py 
and polarization vector e v is similarly parameterized in terms of four form factors V, 
A , Ai and A 2 , again functions of the exchanged momentum squared q 2 = (pn —Pv) 2 - 
Here V denotes the vector form factor and is expected to be dominated by vector 
meson resonance exchange, the axial A\ and A 2 form factors are expected to be 
dominated by axial resonances, while Aq denotes the pseudoscalar form factor and 
is expected to be dominated by pseudoscalar meson resonance exchange [18]. As in 
previous case in order that the matrix elements are finite at q 2 = 0, the form factors 
must also satisfy the well known relation y4o(0) + y4 1 (0)(mH + my)/2my — A 2 (0)(m H — 
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my)/2mv = 0. 

Next we follow the analysis of Ref. [9], where the F + form factor in H — > P 
transitions is given as a sum of two pole contributions, while the Fq form factor is 
written as a single pole. This parametrization includes all known properties of form 
factors at large m#. Using a relation which connects the form factors within large 
energy release approach [19] the authors in Ref. [!9J propose the following form factor 
parametrization 

F+ M = r®— F„tf) = m, (3) 

(1 — x)[l — ax) I — ox 

where x = q 2 /m 2 H *. 

Utilizing the same approach we propose a general parametrization of the heavy 
to light vector form factors, which also takes into account all the known scaling and 
resonance properties of the form factors [151 US] As already mentioned, there exist 
the well known HQET scaling laws in the limit of zero recoil [20J while in the SCET 
limit q 2 — > one obtains that all four H — > V form factors can be related to only two 
universal SCET scaling functions [19J. 

The starting point is the vector form factor V, which is dominated by the pole at 
t = m 2 H * when considering the part of the phase space that is close to the zero recoil. 
For the heavy — > light transitions this situation is expected to be realized near the 
zero recoil where also the HQET scaling applies. On the other hand, in the region 
of large recoils, SCET dictates the scaling described in [19]. In the full analogy with 
the discussion made in Refs. [HI [21], the vector form factor consequently receives 
contributions from two poles and can be written as 

where x = q 2 /m 2 H » ensures, that the form factor is dominated by the physical H* 
pole, while a measures the contribution of higher states which are parametrized by 
another effective pole at m 2 s = m 2 H , / a. 

An interesting and useful feature one gets from the SCET is the relation between 
V and A\ [T9J [221 ESI EI] at q 2 ps 0. When combined with our result (j3J, it imposes a 
single pole structure on A\. We can thus continue in the same line of argument and 
write 

1 — b'x 

Here £ = m 2 H /(mH + my) 2 is the proportionality factor between A\ and V from the 
SCET relation, while b' measures the contribution of resonant states with spin-parity 
assignment 1 + which are parametrized by the effective pole at mL,* = m 2 H ,/b'. It 

eff 

can be readily checked that also Ai, when parametrized in this way, satisfies all the 
scaling constraints. 
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Next we parametrize the Aq form factor, which is completely independent of all 
the others so far as it is dominated by the pseudoscalar pole and is proportional to a 
different universal function in SCET. To satisfy both HQET and SCET scaling laws 
we parametrize it as 



2 A (0) 

(i - - a y) 

where y = q 2 /m 2 H ensures the physical 0~ pole dominance at small recoils and a' 
again parametrizes the contribution of higher pseudoscalar states by an effective pole 
at m 2 H , = m 2 H /a'. The resemblance to V is obvious and due to the same kind of 

eff 

analysis [9] although the parameters appearing in the two form factors are completely 
unrelated. 

Finally for the A 2 form factor, due to the pole behavior of the Ai form factor 
on one hand and different HQET scaling at g^ ax on the other hand, we have to go 
beyond a simple pole formulation. Thus we impose 

, , 2 , A 2 (0) 

Mq } = (l-fs)(l -&»*)' (7) 

which again satisfies all constraints. Due to the relations between the form factors we 
only gain one parameter in this formulation, b" . This however causes the contribution 
of the 1 + resonances to be shared between the two effective poles in this form factor. 

At the end we have parametrized the four H — > V vector form factors in terms 
of the six parameters V(0), A (0), a, a', b' and b" (A 2 (0) is fixed by the kinematical 
constraint). 

In our heavy meson chiral theory (HM%T) calculations we use the leading order 
heavy meson chiral Lagrangian in which we include additional charm meson reso- 
nances. The details of this framework are given in jll] and [15] . We first calculate 
values of the form factors in the small recoil region. The presence of charm meson 
resonances in our Lagrangian affects the values of the form factors at g^, ax and in- 
duces saturation of the second poles in the parameterizations of the F + (q 2 ), V(q 2 ) 
and A (q 2 ) form factors by the next radial excitations of DJ S \ and D( s ) mesons re- 
spectively. Although the D mesons mat not be considered heavy enough, we employ 
these parameterizations with model matching conditions at g^ax- Using HQET pa- 
rameterization of the current matrix elements [EUdS], which is especially suitable for 
HM^T calculations of the form factors near zero recoil, we are able to extract consis- 
tently the contributions of individual resonances from our Lagrangian to the various 
D — ► P and D — > V form factors. We use physical pole masses of excited state 
charmed mesons in the extrapolation, giving for the pole parameters a = m 2 H » / m 2 H „ , 
a' = m 2 1 /m 2 H ,, b' = m 2 H ,/m 2 HA . Although in the general parameterization of the 
form factors the extra poles in F + , V and A),i,2 parametrized all the neglected higher 
resonances beyond the ground state heavy meson spin doublets (0~, l - ), we are here 
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saturating those by a single nearest resonance. The single pole q 2 behavior of the 
Ai(q 2 ) form factor is explained by the presence of a single 1 + state relevant to each 
decay, while in ^(g 2 ) in addition to these states one might also account for their 
next radial excitations. However, due to the lack of data on their presence we assume 
their masses being much higher than the first 1 + states and we neglect their effects, 
setting effectively b" = 0. 

The values of the new model parameters appearing in D — > Plvi decay ampli- 
tudes [H] are determined by fitting the model predictions to known experimental 
values of branching ratios B(D° -> K~l+v), B(D + -»• K°£+is), B(D° -> t~£+u), 
B(D + -> 7T°£ + z/), B(D+ -> t)£+p) and B(D+ -> r/£+z/) [6J. In our calculations of 
decay widths we neglect the lepton mass, so the form factor Fq, which is proportional 
to g M , does not contribute. For the decay width we then use the integral formula 
proposed in [25J with the flavor mixing parametrization of the weak current defined 
in [HI. 

Similarly in the case of D — ► Vlv\ transitions we have to fix additional model 
parameters [15] and we again use known experimental values of branching ratios 
B(D -> K*~£+v), B(D+ -> M + v), B(D + -> p°£ + u), B(D + -> if*°£ + z/), as well as 
partial decay width ratios T L /T T {D+ -> K*°£+is) and r+/r_(£>+ -> K*°£+u) ©. 
We calculate the decay rates for polarized final light vector mesons using helicity 
amplitudes -ff+,-,o as in for example [26]. By neglecting the lepton masses we again 
arrive at the integral expressions from [25J with the flavor mixing parametrization of 
the weak current defined in [T5] . 

We first draw the q 2 dependence of the F + and Fq form factors for the D° — ► 
D° — > 7T _ and D s — > transitions. The results are depicted in Fig. [2j Our 
model results, when extrapolated with the double pole parameterization, agree well 
with previous theoretical [271 EE] and experimental [7J, [8] studies whereas the single 
pole extrapolation does not give satisfactory results. Note that without the scalar 
resonance, one only gets a soft pion contribution to the Fq form factor. This gives for 
the q 2 dependence of F a constant value for all transitions, which largely disagrees 
with lattice QCD results [2SJ as well as heavily violates known form factor relations. 

We also calculate the branching ratios for all the relevant D — > P semileptonic 
decays and compare the predictions of our model with experimental data from PDG. 
The results are summarized in Table [2j For comparison we also include the results 
for the rates obtained with our approach for F + (q 2 aSuX ) but using a single pole fit. It 
is very interesting that our model extrapolated with a double pole gives branching 
ratios for D — > Plv^ in rather good agreement with experimental results for the 
already measured decay rates. It is also obvious that the single pole fit gives the rates 
up to a factor of two larger than the experimental results. Only for decays to r] and 
rf as given in Table an agreement with experiment of the double pole version of 
the model is not better but worse than for the single pole case. 

We next draw the q 2 dependence of all the form factors for the D° — > K~*, D° — > 
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Figure 2: q 2 dependence of the D° — > K (upper left), D° — > n (upper right) and 
D s — > X° (lower) transition form factors. 



Decay 


B (Mod. double pole) 


[%) B (Mod. single pole) 


[%) B (Exp. PDG) [%) 


D° — ► 




3.4 


4.9 


3.43 ±0.14 




7T~ 


0.27 


0.56 


0.36 ±0.06 




7/ 


1.7 


2.5 


2.5 ±0.7 




f/ 


0.61 


0.74 


0.89 ±0.33 


D+ -> 


Z° 


9.4 


12.4 


6.8 ±0.8 


-> 


7T° 


0.33 


0.70 


0.31 ±0.15 


-> 


r/ 


0.10 


0.15 


< 0.5 


-> 


77' 


0.016 


0.019 


< 1.1 




K° 


0.20 


0.32 





Table 2: The branching ratios for the D —> P semileptonic decays. Comparison of 
our model fit with experiment as explained in the text. 

p~ and D s —>■ <p transitions. The results are depicted in Fig. [3j Our extrapolated 
results for the shapes of the D — > V semileptonic form factors agree well with existing 
theoretical studies [23 ESI ED], while currently no experimental determination of 
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Figure 3: q 2 dependence of the D° — > K* (upper left), D° — > p (upper right) and 
D s — > (lower) transition form factors. 

the form factors' shapes in these decays exists. 

We complete our study by calculating branching ratios and partial decay width 
ratios also for all relevant D — > Vtvg decays. They are listed in Table [3] together with 
known experimentally measured values. 



Decay 


B (Mod.) 


[%] B (Exp.) [%] 


r L /r T (Mod.) 


r + /r_ (Mod.) 


D ^K* 


2.2 


2.15 ±0.35 [S] 


1.14 


0.22 


D ^p 


0.20 


0.194 ±0.039 ±0.013 [31] 


1.11 


0.14 


D+ -> K* 


5.6 


5.73 ± 0.35 [6] 


1.13 


0.22 


D+^ Po 


0.25 


0.25 ±0.08 [6] 


1.11 


0.14 


D + -> uj 


0.25 


0.17 ±0.06 ±0.01 [31] 


1.10 


0.14 


D s ^§ 


2.4 


2.0 ±0.5 [6] 


1.08 


0.21 


D s -> K * 


0.22 




1.03 


0.13 



Table 3: The branching ratios and partial decay width ratios for the D — > V semilep- 
tonic decays. Comparison of our model fit with experiment as explained in the text. 
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3 Search for new physics in rare D decays 



At low-energies new physics is usually expected in the down-like quark sector. Nu- 
merous studies of new physics effects were performed in the s —* d, b — > s(d), sd ds, 
bd <-> db and bs <-> sb transitions. 

However, searches for new physics in the up-like quark sector at low energies were 
not so attractive. Reasons are following: a)flavor changing neutral current processes 
at loop level in the standard model suffer from the GIM cancellation leading to very 
small effects in the c — > u transitions. The GIM mechanism acts in many extensions 
of the standard model too, making contributions of new physics insignificant, b) Most 
of the charm meson processes, where c — > u and cu <->■ cu transitions might occur are 
dominated by the standard model long-distance contributions [33] - [41 J. 

On the experimental side there are many studies of rare charm meson decays. The 
first observed rare D meson decay was the radiative weak decay D — ► 7 . Its rate 
BR(D -> 7 ) = 2.6t8;e x !0~ 5 has been measured by Belle collaboration [42J and 
hopefully other radiative weak charm decays will be observed soon|43j. 

In the standard model (SM) [33] the contribution coming from the penguin dia- 
grams in c — > U7 transition gives branching ratio of order 10~ 18 . The QCD corrected 
effective Lagrangian [44] gives BR(c — * wy) ~ 3 x 10~ 8 . A variety of models beyond 
SM were investigated and it was found that the gluino exchange diagrams [45] within 
general minimal supersymmetric SM (MSSM) might lead to the enhancement 

BR(c -» u 7 )mssm _ 1Q 2 ^ 
BR(c -> U7) SM 

Within SM the c — > ul + l~ amplitude is given by the 7 and Z penguin diagrams 
and W box diagram. It is dominated by the light quark contributions in the loop. The 
leading order rate for the inclusive c — > ul + l~ calculated within SM [3S] was found to 
be suppressed by QCD corrections [34] • The inclusion of the renormalization group 
equations for the Wilson coefficients gave an additional significant suppression [40] 
leading to the rates T(c — > ue + e~)/T D o = 2.4 x 10~ 10 and T(c — * ufi + fi~)/T D o = 
0.5 x 10~ 10 . These transitions are largely driven by virtual photon at low dilepton 
mass mil. 

The leading contribution to c —>■ ul + l~ in general MSSM with conserved R parity 
comes from the one-loop diagram with gluino and squarks in the loop [341 13"9l 145] . 
It proceeds via virtual photon and significantly enhances the c — > ul + l~ spectrum at 
small dilepton mass m«. The authors of Ref. [34j have investigated supersymmetric 
(SUSY) extension of the SM with R parity breaking and they found that it can modify 
the rate. Using most resent CLEO [43] results for the D + — > 7i + fi + fi~ one can set 
the bound for the product of the relevant parameters entering the R parity violating 
^22fc^2ifc — 0.001 (assuming that the mass of squark M5 ~ 100 GeV). This bound 
gives the rates BR R (c — > ue + e~) ~ 1.6 x 10~ 8 and BR R (c — ► u[i + [i~) ~ 1.8 x 10~ 8 . 
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Some of models of new physics (NP) contain an extra up-like heavy quark inducing 
flavor changing neutral currents at tree level for the up-quark sector [32 | l4"6l Wf\ |4"9~ 1 150] . 
The isospin component of the weak neutral current is given in [32] as 

*m = \u?^U'l - i5£yz>£ (9) 

with L = ~(1 —75) and mass eigenstates = (uL,CL,t L ,T L ) T , D™ = (d L , s_l, &l) t . 
The neutral current for the down-like quarks is the same as in the SM, while there are 
tree-level flavor changing transitions between up-quarks if O ^ I. The elements of 
4x4 matrix Q can be constrained by CKM unitarity violations currently allowed by 
experimental data. Even more stringent bound on cuZ coupling Q uc comes from the 
present bound on Am in D° — D° transition. It gives \Q UC \ < 0.0004 and we use the 
upper bound to determine the maximal effect on rare D decays in what follows. In 
this case the dilepton mass distribution of the c — > ul + l~ differential branching ratio 
can be enhanced by two orders of magnitude in comparison with SM (see Fig. 4). 

A particular version of the model with tree-level up-quark FCNC transitions is the 
Littlest Higgs model [51]. In this case the magnitude of the relevant c — ► uZ coupling 
Q C u = \V u b\\V c b\v 2 / f 2 < 10~ 5 is even further constrained via the scale / > Oil TeV) 
by the precision electro-weak data. The smallness of Q uc implies that the effect of 
this particular model on c — > ul + l~ decay and relevant rare D decays is insignificant 
[32]. 

The study of exclusive D meson rare decay modes is very difficult due to the 
dominance of the long distance effects [33] - |38| . The inclusive c — > ul + l~ can be 
tested in the rare decays D -> fi + fi~, D -> P(V)l + l~ [Ml EE [35]. 

The branching ratio for the rare decay D — > fx + fi~ is very small in the SM. The 
detailed treatment of this decay rate |34j gives Br(D — > fi + fi~) — 3 x 10~ 13 [31] . 
This decay rate can be enhanced within a study which considers SUSY with R parity 
breaking effects [3U ST]. Using the bound X 22 k^2ik — 0-001 one obtains the limit 
Br(D -> - 4 x !0" 7 - 

The D — > P(V)l + l~ decays offer another possibility to study the c — > ul + l~ 
transition in charm sector. The most appropriate decay modes for the experimental 
searches are D + — ► 7r + / + / _ and D° — > p°e + e~ . In the following we present the possible 
maximal effect on these decays coming from a general class of models with tree level 
cuZ coupling at its upper bound \Q UC \ = 0.0004. We already pointed out that in the 
Littlest Higgs model, which is a particular version of these models, the coupling fl uc 
is constrained to be smaller and the effects on rare D decays are insignificant [32] . 

The calculations of the long distance contributions in the decays D + — > 7r + / + / _ 
and D° — > p°l + l~ are presented in Refs. [321 ESI EH]- The contributions of the 
intermediate vector resonances Vq — P ^^ with Vq — > l + l~ constitute an important 
long-distance contribution to the hadronic decay, which may shadow interesting short- 
distance contribution induced by c — > ul + l~ transition. 
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Figure 4: The dilepton mass distribution dBr /dm 2 ee for the inclusive decay c — > ul + l 
as a function of the dilepton mass square m 2 ee = (p + +P-) 2 . 

Our determination of short and long distance contributions to D + — ► 7r + / + /~ takes 
advantage of the available experimental data [32j. This is a fortunate circumstance 
for this particular decay since the analogous experimental input is not available for 
determination of the other D — > Xl + l~ rates in a similar way. The rate resulting 
from the amplitudes (14) and (19) of [32] with \VL UC \ = 0.0004 are given in Fig. 5 and 
Table 4. 

We are unable to determine the amplitude of the long-distance contribution to 
D° — > p°Vo — > p°l + l~ using the measured rates for D° — ► p°Vo since only the rate of 
D° — > p°(f) is known experimentally. We are forced to use a model [38], developed to 
describe all D — > V^Z + /~ and D — > V7 decays, and the resulting rates are presented 
in Fig. 6 and Table 4. 

Therefore, the total rates for D — ► Xl + l~ are dominated by the long distance 
resonant contributions at dilepton mass m« = m p , ra w , and even the largest 
contributions from new physics are not expected to affect the total rate significantly 
[341 [39] . New physics could only modify the SM differential spectrum at low mu 
below p or spectrum at high mu above 0. In the case of D — ► ttI + 1~ differential decay 
distribution there is a broad region at high mu (see Fig. 5), which presents a unique 
possibility to study c — ► ul + l~ transition [39, 132]. 

The non-zero forward-backward asymmetry in D —>■ pl + l~ decay arises only when 
C10 7^ (assuming mi — > 0). The enhancement of the C10 in the NP models [32] is 
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new physics: tree-level c L ->u L Z coupling (O =0.0004) 
+ + + - 

D -> 71 e e 

10- 3 , 




12 3 



m 2 [GeV 2 ] 



Figure 5: The dilepton mass distribution dBr /dm 2 ee for D + — > n + e + e . 

due to the tree-level u^^c^Z^ coupling and leads to nonzero asymmetry ApB^m^) 
shown in Fig. 7. The forward-backward asymmetry for D° — > p°l + l~ vanishes in SM 
(Cio — 0), while it is reaching O(10~ 2 ) in NP model with the extra up-like quark as 
shown in Fig. 7. Such asymmetry is still small and difficult to be seen in the present 
or planned experiments given that the rate itself is already small. 

We have investigated impact of the tree-level flavor changing neutral transition 
c — > uZ on the rare D meson decay observables. However, the most suitable D + — > 
7T + / + /~ and D° — > p°l + l~ decays are found to be dominated by the SM long distance 
contributions. Only small enhancement of the differential mass distribution can be 
seen in the case of D + — > decay at high dilepton mass and tiny forward 





Br 




short distance 


total rate ~ 


experiment 








contribution only 


long distance contr. 




SM SM + NP 


D+ - 






6 x KT 12 8 x l(T y 


1.9 x 10^ 


< 7.4 x 1(T S 


D+ - 


-> 7T + fl 


V" 


6 x 10" 12 8 x 1(T 9 


1.9 x 10" 6 


< 8.8 x 1(T 6 


D°- 




e" 


negligible 5 x 10 -10 


1.6 x 1(T 7 


< 1.0 x 10" 4 


D° - 




p" 


negligible 5 x 10 -10 


1.5 x 1Q- 7 


< 2.2 x 1Q- 5 



Table 4: Branching ratios for the decays in which c — > ul + l transition can be probed. 
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Figure 6: The dilepton mass distribution for D° p°e + e 
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Figure 7: The forward-backward asymmetry for D° — > p°e + e . 
backward asymmetry can be induced by new physics in D° — > p°l + l~ decay. 
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We conclude that the NP scenarios which contain an extra singlet heavy up-like 
quark, have rather small effects on the charm meson observables. 
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